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Abstract 

Extending previous QCD Hamiltonian studies, we present a new renormaliza- 
tion procedure which generates an effective Hamiltonian for the gluon sector. 
The formulation is in the Coulomb gauge where the QCD Hamiltonian is 
renormalizable and the Gribov problem can be resolved. We utilize elements 
of the Glazek and Wilson regularization method but now introduce a continu- 
ous cut-off procedure which eliminates non-local counterterms. The effective 
Hamiltonian is then derived to second order in the strong coupling constant. 
The resulting renormalized Hamiltonian provides a realistic starting point for 
approximate many-body calculations of hadronic properties for systems with 
explicit gluon degrees of freedom. 
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I. INTRODUCTION 



A quantum field theory that is local will ensure that causality is not violated, however, 
application of such a theory in four- dimensional space-time will generate divergences. About 
50 years ago in a series of papers |T| , Dyson developed renormalization to specifically address 
such infinites in QED. He ingeniously converted the divergence problem into a useful crite- 
rion for selecting appropriate theories, namely that they should be renormalizable. Today 
the view towards renormalization is somewhat more practical and perhaps limited in scope. 
Guided by Wilson and Weinberg, theories are now regarded more as partial descriptions 
rather than complete constructs. Accordingly, contemporary approaches embody the con- 
cept of an effective theory that only includes the necessary degrees of freedom appropriate 
to a specific energy range. 

In this paper, we follow this philosophy and consider the development of an effective QCD 
Hamiltonian for systems with explicit gluon excitations in the few GeV energy range. This 
work extends our previous Hamiltonian glueball investigation |2| and recent renormalization 
treatment of the quark sector [[3]. Our motivation is three-fold. First, we seek a rigorous 
formalism connecting QCD to effective and tractable model Hamiltonians. Second, we wish 
to develop a realistic theoretical treatment to accurately describe hadron structure, especially 
systems with gluonic degrees of freedom such as glueballs and hybrids. Finally, we want to 
understand the basic confinement and chiral symmetry breaking mechanisms. 

We begin with the exact QCD Hamiltonian formulated in the Coulomb gauge and fo- 
cus upon the pure gluonic sector. Renormalizing the quark Hamiltonian, which has been 
previously studied using a sharp cut-off regularization ||, and the quark-glue sector will 
be addressed in a future publication. There are several reasons for utilizing the Coulomb 
gauge in which the divergence of the color vector potential vanishes (V ■ A = 0). As detailed 
by Zwanziger [[|, not only is the Hamiltonian renormalizable in this gauge but the Gribov 
problem || can also be resolved. The essence of the Gribov problem is that specification 
of V ■ A = does not uniquely fix the gauge in non-Abelian gauge theories. In general, 
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there are many copies of each gauge field configuration, all with the same divergence, that 
are related by gauge transformations. As the true physical configuration space of a gauge 
theory is the set of gauge potentials modulo local gauge transformations, one must select a 
single representative from each set of gauge-equivalent configurations. The resulting sub-set, 
-4phys, of independent field configurations is known as the fundamental modular region. 

A convenient characterization of *4. p h ys is given by the "minimal" Coulomb gauge, ob- 
tained by minimizing a suitably-chosen functional along gauge orbits. We define the L2 
norm of a field configuration along a gauge orbit by 



Then *4. p h ys is specified by choosing from each gauge orbit the configuration which globally 
minimizes F. It is straightforward to verify that the extrema of F satisfy V • A = 0, so 
that this defines a Coulomb gauge. Furthermore, at a minimum of F the Faddeev-Popov 
operator satisfies M = — V-D > (i.e., its eigenvalues are non-negative). Here D = V — igA 
is the covariant derivative. One can show that configurations for which M = occur only 
on the boundary of -4 p h ys - These actually represent gauge copies and must be identified to 
give the full physical configuration space. Needless to say, the resulting functional space has 
an extremely complicated structure. Zwanziger has shown how to rigorously implement the 
restriction to *A p h ys in the infinite-volume limit for a gauge-fixed version of the lattice QCD 
Hamiltonian formulation [|J. 

The confinement phenomenon in QCD has two complementary aspects: (1) there is a 
long range attractive potential between colored sources; (2) the gluons which mediate this 
force are absent from the spectrum of physical states. This poses something of a mystery 
in a covariant gauge, however, since for small q 2 , (1) suggests that the gluon propagator 
should be more singular than 1/q 2 , while (2) implies the propagator is suppressed. Thus the 
mechanism for confinement is not particularly transparent [0]. 





where G(x) is a gauge transformation and 



A G = GAG 1 " - GVG f . 
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In the Coulomb gauge, in contrast, these two aspects can comfortably co-exist: the long 
range force is represented by the instantaneous Coulomb interaction, while the (transverse) 
gluon propagator can simultaneously be suppressed at q 2 = 0. Indeed, a detailed analysis 
suggests precisely this picture. The key issue is the proper identification of the fundamental 
modular region. Zwanziger has argued recently that restricting the configuration space to 
the fundamental region results in the Coulomb term acquiring a singular contribution at 
long distances |§ .[] This singularity is also connected to a suppression of the propagator for 
the would-be-physical transverse gluons at q 2 = 0. This suggests a rather appealing scenario 
for confinement in the Coulomb gauge. 

For these reasons and because we are interested in understanding confinement, we adopt 
the Coulomb gauge and only work with transverse components of the gluon fields. Next 
we divide the Hamiltonian H into a free part, H Q , which is the full Hamiltonian evaluated 
with zero coupling constant, and interaction defined by Hi = H — Hq. We work in a Fock 
space spanned by eigenfunctions of H Q with eigenvalues E n . In general matrix elements of 
the interaction diverge in this space and we regulate by suppressing contributions between 
states in which the energy difference \E n — E m \ is larger than scale set by a cut-off pa- 
rameter A. However, we do not use the sharp cut-off procedure of Glazek and Wilson |3j 
since this generates unacceptable non-local interactions upon renormalization. These may 
be avoided by using a smooth regulator. This interesting result is not attributed to the 
specific Hamiltonian or choice of gauge since this feature emerges in scalar field theories as 
well. Although our theory is now rendered finite it is inappropriately dependent upon the 
parameter A. Furthermore, since this parameter is governed by free energies, the theory 
is no longer Lorentz invariant. It is also not gauge invariant but after renormalizing both 
symmetries will be restored along with elimination of all the cut-off dependence. In addi- 
tion to this regularization we also suppress divergent matrix elements of one-body operators 



This mechanism for confinement was first suggested by Gribov ||. 
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which arise from normal ordering the two-body interaction. Again a continuous, exponential 
cut-off regulator is used and cut-off sensitivity is removed by renormalization. 

Renormalization is achieved by adding a counterterm Hamiltonian which can be ex- 
pressed in terms of unknown coefficients and a complete set of local operators that respect 
the symmetries of the regularized Hamiltonian. An effective and elegant means for deter- 
mining these coefficients is by performing similarity transformation ||. Because the trans- 
formation is unitary, the physical content is preserved but the cut-off is now reduced to a 
lower value Ai. The transformation also introduces new interactions which incorporate the 
physics contained between the scales A and Ai. A key aspect of the renormalization scheme 
is the requirement that the transformed Hamiltonian be form invariant, i.e. it maintains its 
mathematical structure but with Ai now replacing A ||. 

Because A can be arbitrarily large, we select an initial scale which is amenable to pertur- 
bation theory. The regularized Hamiltonian is expanded in powers of the coupling constant 
g and for this study only terms to order g 2 are retained. While the similarity transformation 
rigorously evolves the scale it cannot account for confinement at hadronic energies. This 
constitutes our main theoretical omission. For hadronic applications confinement may be 
described by supplementing the renormalized Hamiltonian with a confining interaction. We 
will address this issue further in a future publication. 

At this point we have obtained an effective Hamiltonian suitable for application to the 
vacuum and excited gluonic states, i.e. glueballs. However, our experience indicates 
that an improved description can be obtained by performing an additional similarity trans- 
formation to a quasiparticle basis. This is a BCS (Bogoliubov) rotation which mixes the bare 
parton gluon creation and annihilation operators. As previously demonstrated |2||| through 
variational calculations for an unrenormalized Hamiltonian in both the quark and gluon 
sectors, dynamical chiral symmetry breaking occurs and reasonable values for condensates 
and constituent masses are obtained. Similar vacuum BCS variational calculations using our 
renormalized effective gluon Hamiltonian are in progress and will be reported in a future 
publication along with many-body Tamm-Dancoff and Random Phase Approximations for 
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the excited hadron (glueball) spectrum. 

This paper is divided into five sections and one appendix. In the next section our nota- 
tion is established and the canonical QCD Hamiltonian is specified. Section III addresses our 
regularization scheme. Section IV presents the main result of the paper. It includes a discus- 
sion of the counterterm Hamiltonian and the form invariance renormalization scheme. This 
is followed by a sub-section devoted to the similarity transformation which yields the final 
renormalized effective interaction. A detailed representation of the canonical Hamiltonian 
in Fock space is presented in the appendix. 



Our starting point is the canonical QCD Hamiltonian in the Coulomb gauge V ■ A = 
[[L(J. We denote spatial vectors by bold-faced quantities and use the matrix notation 
for the gauge fields: A = A a T a , with T a the generators of SU(iV c ) in the fundamental 
representation, satisfying [T a ,T b ] = if abc T c . The gauge covariant derivative is D = V — igA 
operating on objects in the fundamental representation and D = V — ig[A, ] for objects in 
the adjoint representation. 

The dynamical degrees of freedom are the transverse gauge fields A, their conjugate 
momenta II (also transverse) and the quark field if). The canonical Hamiltonian takes the 
form 



II. CANONICAL HAMILTONIAN 



H, 



can 



H q + H g + H c + H t 



(2.1) 



with 




(2.2) 




(2.3) 




(2.4) 
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H qg — —g J d 3 x^(-x)ct ■ A(x)-0(x). 



(2.5) 



Here a 1 = 7°7* and f3 = 7 are the Dirac matrices, J is the Faddeev-Popov determinant, 



J = det [V • D] 



(2.6) 



normalized so that det [V 2 ] = 1, and B is the non-Abelian magnetic field with components 



B? = eqkVjAl + 9 -e l]k f abc A)Al 



(2.7) 



In the Coulomb term, Eq. fl2.4|) , the kernel K is represented in "matrix" notation as 
K o6 (x, y) = (x, a|(V • D)- 1 (-V 2 )(V • D)" 1 ^, b) 



(2.8) 



and the color charge density p a includes both quark and gluonic contributions 



p a (x) = ^ (x ) T>(x ) + /^A b (x) • n c (x). 



abc \ bi 



(2.9) 



The nonvanishing canonical commutator for the gauge field is 



ViV, 



A»(x), n*(y) = tS ab ( By - ^ ) 5 (3) (x - y) 



(2.10) 



The Fourier decompositions are 

d 3 k 1 



A«(x) 



(2vr) 3 



<(k) + af(-k) 



,t'k-x 



(2.11) 



n?(x) = -i I ^ 



(2.)3V 2 l« k )-«? , (- k ) 
with c<Jk = |k| = k. In momentum space Eq. (|2.10| ) reduces to 



_ikx 



(2.12) 



<(k),af (k')l = 5 ab (2vr) 3 5( 3 )(k- kOAi(k), 



(2.13) 



where 



A,(k)^-^. 



(2.14) 



Note that transversality of A and II requires 



k ■ a a (k) = k ■ a at (k) = 0. (2.15) 
These operators can also be decomposed in terms of orthogonal polarization vectors 



a? 



(k)= J2 Q(k,AK(k,A). (2.16) 



A=l,2 



For the Fermi field the operator expansion takes the form 

<M X ) = E / 7^ H k > *) fo (k, *) + v(-k, s)J(-k, s)] e ik ' x , (2.17) 

where s labels the helicity (color and flavor indices have been suppressed) and a nonrela- 
tivistic normalization is used such that v)u = v'v = 1. The nonvanishing canonical anti- 
commutators are then 

{6(k, s), &+(k', s')} = {d(k, s), s')} = 5 s A2n) 3 5^(k - k'). (2.18) 

It is straightforward to express the canonical Hamiltonian in the Fock representation. 
The result to 0(g 2 ) is given in the Appendix. 



III. REGULARIZATION SCHEME 

Perturbative schemes for renormalizing Hamiltonians typically suffer from technical prob- 
lems related to the occurrence of vanishing energy denominators. An elegant way of avoiding 
this difficulty is the cut-off method of Glazek and Wilson . This approach uses the basis 
formed by eigenstates of the free Hamiltonian H , with eigenvalues E n . The theory is then 
regulated by suppressing matrix elements of H between states for which the difference in free 
energies E nm = E n — E m is large. Specifically, matrix elements of the regulated Hamiltonian 
H(A) are defined by 

(n\H(A)\m) = E n 6 nm + / nm (A)(n|F 7 |m), (3.1) 
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where Hi = H — Hq is the interaction. Here f nm (A) is some convenient function which 
approaches unity for \E nm \ <C A and vanishes for \E nm \ ^> A. In this work we employ a 
smooth cutoff of the formP] 

/ nm (A)=e" E -/ A2 . (3.2) 

Note that diagonal matrix elements are unmodified. 

A few remarks on this regularization scheme are in order. First, it does not reduce the 
size of the Hilbert space, that is, no states are removed from the theory. Reduction to a finite 
number of degrees of freedom — which is necessary for the nonperturbative calculations we 
wish eventually to perform — occurs later, after the variational vacuum state is determined 
and the transition to constituent quarks and gluons has been achieved. At this point we can 
apply, for example, a Tamm-Dancoff truncation and consider hadronic states built from a 
small number of constituents. This truncation can perhaps be justified by the emergence of 
constituent- scale masses for the effective degrees of freedom, which kinematically suppresses 
mixing with more complex, multi-constituent states. A similar truncation in the bare parton 
basis (current quarks and mass zero gluons) would be essentially meaningless. 

Second, since the cut-off is defined in terms of free energy differences it is not fully Lorentz 
invariant (though it is rotationally invariant, of course). This is unavoidable in a Hamiltonian 
framework, particularly if positivity of the Hilbert space norm is to be preserved (thus ruling 
out, e.g., the Pauli-Villars scheme). Indeed we conjecture that requiring positivity forces 
one to employ a noncovariant gauge. The renormalized effective Hamiltonian will therefore 
contain Lorentz-noninvariant operators, that is, operators which correspond to Lorentz- 
noninvariant terms at the Lagrangian level. The regulator also violates gauge invariance, 



2 We have also investigated a sharp cut-off, with f nm = 6 (A — \E nm \). However, this choice leads 
to pathologies in the renormalized Hamiltonian, specifically nonlocal counterterms. These terms 
also arise in scalar field theory and thus have nothing to do with the choice of the Coulomb gauge. 
This issue is discussed further below. 
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though this may be of less significance since we work in a fixed gauge. Nevertheless our 
renormalization procedure, to be detailed below, will automatically provide the counterterms 
necessary to restore both symmetries. 

Finally, this scheme does not completely regulate the theory. Normal ordering the two- 
body operators in the canonical Hamiltonian leads to one-body operators (see Eqs. ( |A4j ) and 
( |A8| ) in the appendix) with divergent matrix elements that are not regulated by the above 
procedure. For these matrix elements we insert an additional cut-off factor, exp(— 2o;^//A 2 ), 
in the integrands of Eqs. flA4]) and (|A8|) -PI While there is substantial freedom in implement- 
ing this supplementary regularization, observables computed from the Hamiltonian will be 
independent of the cut-off prescription. We remark on this below. 

At this point we have a fully regulated Hamiltonian formulation of QCD in the Coulomb 
gauge. The next task is to remove the dependence on the cut-off parameter A by adding 
counterterms to the Hamiltonian. This will be done perturbatively, which should be reason- 
able for QCD if the cut-off is not too low. The resulting renormalized effective Hamiltonian, 
-ffeff(A), can then be analyzed nonperturbatively using many-body techniques. 

IV. CUT-OFF DEPENDENCE AND RENORMALIZATION 

A. Counterterm Hamiltonian and Form Invariance 

Our goal is to construct the renormalized effective Hamiltonian, H e g(A), which contains 
a cut-off but gives cut-off independent results through some fixed order in perturbation 
theory. It can be expressed as the sum of the canonical Hamiltonian and a "counterterm" 
Hamiltonian 



Of course, we could simply omit these operators from the Hamiltonian since our goal is to remove 
the cut-off dependence by adding counterterms. However, the divergences from these terms cancel 
divergences elsewhere in the theory so we retain them. 
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ifeff(A) = ^can(A) + #Ct(A). (4.1) 

The canonical Hamiltonian contains normal ordering contributions with integrals regulated 
as previously described. The counterterm Hamiltonian, which begins at order g 2 in pertur- 
bation theory, has operators with the same structure along with additional terms necessary 
to correct for violations of Lorentz invariance introduced by the regulator. Note also that 
the matrix elements of Hqt depend explicitly on A in addition to containing the regulating 
exponentials. In general, 

tf CT (A) = ]>XA)A^ A , (4.2) 

i 

where the {Of} are a complete set of local operators invariant under the symmetries pre- 
served by the regulator. The superscript A indicates that the regulating functions are to be 
associated with the operators Oi] that is, 

(n\Of\m) =e- E ^ /A2 (n\Oi\m), (4.3) 

where Oi is the "bare" operator. The explicit powers of A are inserted in Eq. (|4.2|) so 
that the coefficients q(A) are dimensionless (jii is four minus the mass dimension of the 
operator Oi). Operators for which n« < are irrelevant in the renormalization group sense 
and will generally be ignored. Note that the coefficients q can depend on A only through 
their dependence on the "canonical" couplings of the theory. Thus 

q(A)=^ 2 (A) + ---, (4.4) 

where the ellipsis indicates higher order terms in the perturbative expansion and for sim- 
plicity only a single independent coupling has been assumed. The pure numbers, q, are 
independent of A. Thus the coefficients depend logarithmically on A and one may only 
specify the change in the coefficient. This may be contrasted with, for example, a gluon 
mass counterterm which arises solely because the regulator violates gauge invariance. In 
this case the gluon mass is completely calculable within perturbation theory. 
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An exhaustive way of determining Hqt would be to list all permissible operators and then 
extract the Cj by requiring observables calculated perturbatively be both cut-off independent 
and Lorentz covariant. A more appealing, physically equivalent, approach adopted here is 
to perform a similarity transformation on the Hamiltonian which reduces A to a new cut- 
off Ai in the regulating exponentials. Because the transformed Hamiltonian is equivalent 
to the original one, calculated observables (e.g., eigenvalues) are unchanged. In addition 
to replacing A — > Ai in H caii) the transformation also generates new, effective interactions 
representing the physics between the scales A and Ai. Schematically, we have 

# can (A) - #can(A!) + SH{A U A) (4.5) 

under the similarity transformation. The counterterm Hamiltonian is then determined by 
requiring that the full, transformed Hamiltonian have the same form as the original one, but 
with A everywhere replaced by Ai. This requirement of form invariance, or coherence, of 
the Hamiltonian under changes of the cut-off ensures that we have the most general Hamil- 
tonian allowed [§,[[]]]• In the language of the renormalization group (RG), the self-similar 
Hamiltonian exists on an infinitely long RG trajectory and thus constitutes a renormalized 
Hamiltonian. 

We now sketch this procedure to first order in a s = g 2 /4:Tr. In this case Hqt is of order 
a s and since the similarity transformation leaves Hqt unchanged to this order (apart from 
replacing A by Ai in just the regulating functions) we have 

H cS (A) - # can (Ai) + SH(A 1 , A) + F CT (A). (4.6) 

In this expression Hqt is understood to have A replaced by Ai in the regulating functions 
but nowhere else. We now demand that the transformed Hamiltonian be equal to 

Fe^Ax) = H cm (A 1 ) + H CT (Ai). (4.7) 

This uniquely determines the structure of Hqt- 

A prescription for treating the cut-off dependent terms from normal ordering must also 
be specified since these are not regulated by the cut-off on energy differences and are also 
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unaffected by the similarity transformation. Our prescription is that upon performing the 
similarity transformation, we replace A by Ai in these terms in the canonical Hamiltonian 
and include the difference between this and the original terms — the "slice" coming from 
momenta between Ai and A — in SH. 

We first determine contributions to 5H coming from the normal-ordering terms. The 
one-body operators in the effective Hamiltonian have the general form 

^ )(A) = / { [ W *M** + n ?( k )] < f (k)a?( k ) 

+iW£ 6 (k)e-( 2 ^ 2 /A 2 [ a ?(k)o5(-k) + h.c] } . (4.8) 

Note the presence of the regulating function in the second term. From rotational invariance, 

II#(k) = A ab (k)5 i3 + B ab (k)k i k j , (4.9) 

with similar structure for . Due to transversality, however, only the term proportional 
to Sij need be considered. From Eqs. ( |A^) and ( |A8|) we can identify the contributions to 
and M?j arising from normal ordering the two-body operators. From the Coulomb term 
we have 

nf a \k) = ^N c 5 ab I ^V(k + k') ( L d±^K\ Dij (k')e-H^ 2 (4.10) 
J 4 J (2ny \ uj k uj k > J 

and 

Mf a \k) = ^N c S« b [ p-V(k + kO ( D^)e-H^. (4.11) 
Next the cut-off is reduced by an infinitesimal amount, A 2 — > A 2 , with 

A 2 = (l-e)A 2 ; e > 0. (4.12) 

The difference 

5U^ )ah (k] A, e) = ng 7 )ab (k; A) - U$p a \k; A x ) (4.13) 



13 



to be included in SH is then expanded to first order in e, followed by an expansion in powers 
of k/A. Retaining only those terms which correspond to relevant or marginal operators in 
the renormalization group sense yields Q 



—NSS ab 
4tt c « 



e 



-A 2 H k 2 + ■ 

3 15 



(4.14) 



where the dots represent terms of order k/A and higher. Similarly we obtain 



5M^ )ab (k-A,e) = ^NA,6 ab { — 



a. 



ab 



An 



(4.15) 



Performing the same steps for the divergent contributions in H g [Eq. ( |A4| )1 produces 



o:. s 



6U^ ao (k;A,e) = ^N c 5 l ,5 



ab 



-A 2 + ■ ■ 



(4.16) 



and 



5Mlf a \k-A 1 e) = ^NA,5 ab [-\r-A 



a 



4:71 



L0 k 



(4.17) 



B. Similarity Transformation 

We now develop the similarity transformation renormalization procedure. This follows 
the formulation of Wegner |12|] who derived flow equations for diagonalizing a Hamiltonian. 
Consider a continuous transformation governed by a unitary operator U which depends 
upon a (flow) parameter A. Let T(A) be the generator of the corresponding infinitesimal 
transformation. Then for a continuous transformation from Aj to Xf 

C/(A /; A,)=e/>:' T<W \ (4.18) 

where A-ordering of the exponential is understood. Note that T must be anti-hermitian 
since W = U -1 . The similarity transformation for a Hamiltonian which also depends on A 
is then given by 



4 It may be necessary to retain some irrelevant operators if the cut-off is reduced much below the 
scale of physical interest, but for the present we ignore these. 
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H(X f ) = U(X f ;X t )H(X t )U\X f ;X i 



(4.19) 



and for an infinitesimal evolution 



tiff (A) 
dX 



[T(X),H(X)]. 



(4.20) 



Wegner has determined that the generator choice 



T(X) = [H ,H(X)], 



(4.21) 



suppresses the off-diagonal matrix elements thereby rendering H "more diagonal" as A in- 
creases. We have independently obtained the same generator choice by requiring that the 
transformation properly change the cut-off parameter in the regulating functions for all ma- 
trix elements of our Hamiltonian. We now demonstrate this but first note, from Eqs. ( 4.20| , 
4.21), the dimensions of A are [E~ 2 ] so to work directly with our energy cut-off parameter 
we substitute A — > A~ 2 . Then combining the above equations yields 

dH{A) 



dA~ 2 



[H ,H(A)],H(A) 



(4.22) 



By expanding the interaction Hamiltonian in powers of the coupling at this scale A, 



ff,(A) = EV 

p=l 



(4.23) 



we can examine Eq. ( 4.22 ) order by order. At zeroth order 

dH 



dA- 



0. 



(4.24) 



so that the transformation preserves the cut-off independence of the free Hamiltonian. At 
first order we find 



rffei(A) 
dA- 2 



[Ho, hi], H 



(4.25) 



Expressed in terms of matrix elements this becomes 

d 



dA~ 2 



(n\h x (A)]m) = -^ m (n|/n(A)|m). 



(4.26) 
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Integrating Eq. ( |4.26| ) from A to Ai then yields 

(nl/ixCAOlm) = e- s -( 1/A ?- 1/A2 )(n|/ il (A)|m). (4.27) 

Since 

(n|/ii(A)|m) cxe~ E ^ A2 , (4.28) 

we see the the transformation has indeed replaced A — > Ai in the regulating exponentials 
for these terms. 

The expression at second order reads 

-jL(n\h 2 {A)\m) = -E 2 nm (n\h 2 (A)\m) 

+ J2(E n + E m - 2E l )(n\h l (A)\l)(l\h l (A)\m). (4.29) 
l 

Once again the first term on the right-hand side has the effect of replacing A — > Ai in 
the overall regulating exponential. It is trivial to see that this same structure occurs for 
each order, so that the transformation does indeed replace A — > Ai in all matrix elements 
of the Hamiltonian. It is convenient to isolate this explicit cut-off dependence by defining 
"reduced" matrix elements, denoted by V p 

(n\h p (A)\m) = e~ E ™ /A * (n\V p {A)\m). (4.30) 

Notice that V\ does not depend on A. In terms of V 2 , Eq. ( |4.29|) reads 



-^ m /A 2 ^L (n |y 2 (A)| m ) = J2(E nl + E^e-^^^inlV^mim). (4.31) 



dA 



-2 



1 

This can be integrated to give 

(n\8V 2 \m) =Y.®(Eni,E ml )(n\V 1 (A)\l}(l\V 1 (A)\m), (4.32) 
1 

where 

(n\8V 2 \m) = (n|y 2 (Ai)|m) - (n|V^(A)|m). (4.33) 
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This change in the interaction explicitly represents the physics removed in reducing the 
cut-off from A to Ax, to this order in perturbation theory. The quantity 

e(E nl , E ml ) = ± (J_ + J_) (e- 2E «^/ A2 - e - 2 ^-/ A ?) (4.34) 

plays the role of an energy denominator in conventional perturbation theory. Note that 
cancellation between the two exponentials suppresses the contribution when the energy 
differences are small. Thus the potentially singular factors 1/E n i, etc., are not problematic. 
It is straightforward to extend this procedure to higher orders. 

We can now compute the various one-body operators induced by the similarity trans- 
formation. Consider first the contribution to 511 arising from qq intermediate states. It is 
given by 

v Jf (-q,s')a j u(k + q,s) , 

"~ K s,s' " V*"' > 

(4.35) 

where 

A = uj k - E q - E k+q . (4.36) 

We again expand this for a small change in the cut-off and discard contributions correspond- 
ing to irrelevant operators. The result is 

5u i,uar k ,i )ab{k . Aj e) = -*L Nf5ij fi>* (£) (1 A 2 + ^kA + .-.), (4.37) 

where the dots represent terms of order 1/A and Nf is the number of light (i.e., lighter than 
the cut-off scale) quark flavors. There is also a contribution arising from qqgg intermediate 
states (analogous to a "Z-graph" contribution) which gives 

su i,uar k ,2 )ab{k . Aj £) = _ZL NfSijl y* (_£_) (_L A 2 ~ ^kA + ■ ■ ■) . (4.38) 
The complete contribution arising from multi-quark intermediate states is thus 
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5U 



(quark,l)ab 
ij 



(k) 



u q 



The contribution from multi-gluon intermediate states may be evaluated similarly; the 
result is 



O'.s 



5U^ lue)ab (k;A,e) = -^N c 5 t ,5 a 



1 A2 !6 72 

, , —A 2 H A; + 

V12 15 



(4.40) 



Combining these results with the contributions from the normal-ordering terms gives the 
complete 511 at this order: 

2Nt 



flI#(*;A,e) 



167T 



H_z^Z A 2 



(4.41) 



where we have set N c = 3. We recognize the lowest-order coefficient of the QCD beta 
function fa = 11 - 2N//3.Q 



Next we consider the various contributions to 5M? . The qqg intermediate states yields 

.uT k J V 12 



<SMt ark)a \k] A, e) = -^-NfSijS^ (—) (^A 2 + •••). (1.12) 



while three-gluon intermediate states give 



,u k J V24 15 

Combining these results with the contributions from the normal-ordering terms gives 

1 



(4.43) 



SM^(k ] A,e) = ^-S tJ 5 ah 



e 



-/3 A 2 - 8r + 



(4.44) 



for N c = 3. 

One-body operators in the counterterm Hamiltonian have the general form 
#ct(A) = J -j^L (u CT (k; A)af(k)<(k) + M CT (k; A)e~^ 2 / A2 [<(k)<(-k) + h.c.]) , 



(4.45) 



since the cut-off dependent terms are all proportional to 5ij5 ab . Let us focus first on the 
term in Hct that depends quadratically on A. On dimensional grounds one has 



5 This is purely accidental since an alternative regulator choice in the normal-ordering terms gives 
a different result. 
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where b is a numerical constant independent of A. After the similarity transformation we 
have 



# can (A) + if CT (A) - H c ^i) + 6H(A l , A) + H CT (A), (4.47) 

where 5H represents the contributions computed above. Note that to this order, the only 
change to if can (A) is to replace A with Ai. Equating the second and third terms on the RHS 
of Eq. (pU7l) to #ct(Ai) yields, using Eq. (g4| ), 

\ef3 A 2 + bA 2 = bAl 



6A 2 (l-e). (4.48) 



Thus b = —/3 /2, and 



ffcT(A) ._-M!/|l(_L) af(kK(k) + .... ( , 49) 

It may be verified that this counterterm removes the quadratically divergent part of the 
gluon self energy to second order, as it should. Of course, this is simply (part of) a gluon 
mass term with the mass given by 

2 «sA)A 2 . . 

m g = g^—- ( 4 - 5 °) 

Note that the quadratically divergent part of SM has exactly this structure. Therefore, 
when expressed in terms of the field operators the complete quadratically divergent part of 
the counterterm Hamiltonian is simply 

2 

#ct(A) = -f j d 3 xA?(x)At(x) + -... (4.51) 

The remaining operators in i^cT may be similarly constructed. Dimensionally, the term 
in Eq. (|4.44D containing the factor —8k 2 represents a logarithmic divergence. Hence, it 
should correspond to the only logarithmic divergence to which we are sensitive at this order, 
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namely the gluon wavefunction renormalization. To see this we invert the Fourier expansions 
and express this operator in terms of the fields to obtain 

a s [AC A 



5H 



Tr J d 3 x (B 2 - n 2 ) , (4.52) 



4tt V 3 

where we have restored the dependence on the number of colors (Ca = N c is the Casimir 
invariant of the adjoint representation of SU(iV c )). This has precisely the correct structure 

1/2 

for a wavefunction renormalization since, if the gauge field is rescaled by A — > Z A A then 
the conjugate momentum is rescaled by the inverse factor, II -> Z A 1/2 Tl. 
The counterterm Hamiltonian will thus contain a term 

H CT = c(A)Tr J d 3 x (b 2 - n 2 ) | s=0 (4.53) 

at this order, with c(A) logarithmically dependent on A. The requirement of coherence of 
the Hamiltonian under a change of the cut-off from A — > Ai then gives 

-^( 1 f i )+ C (A) = c(A 1 ). (4.54) 

Expanding c(A x ) about A x = A yields 

A-^c(A) = ^ ( 8 -^) (4.55) 



dA y ' 4tt V 3 

to 0(e) and hence 



c(A) = M^WA/Ao), (4.56) 



4tt V 3 

where Ao is an arbitrary scale. 

As discussed above, one can interpret this counterterm as corresponding to a renormal- 
ization of the gauge field with form 



A' 



1 + ~c(A) 



A = Z A /2 A. (4.57) 



3 A quick way to see this is to note that in a functional Schrodinger representation II = — ibjbA. 
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Thus, the free part of H can combines with this part of the counterterm Hamiltonian to 
simply produce a free Hamiltonian when written in terms of A'. The gauge field anomalous 
dimension is now easily computed to be 



1a 



--A^-lnZ A 
2 dA 



9 2 (4C. 



16tt 2 



(4.58) 



At lowest order, there is no contribution to 7^ from fermions in this scheme. 

It is significant to note that with a sharp cut-off, the above results are again obtained, 
but with two differences. One is a slight difference in the value of the gluon mass. The 
other, which is problematic, is a contribution 

d 3 k 



Hnl=[ ^) £PoA J (2tt)3 



ai(k)a,i(-k) + o|(k)a|(— k) 



(4.59) 



to the counterterm Hamiltonian. This corresponds to a nonlocal interaction, as may easily be 
seen by inverting the field expansions and writing Eq. ( 4.59 ) in terms of the field operators 

1 



H 



NL 



47T 



e/? A / d 3 xd 3 y 



l x -y| 



A(x) ■ V 2 A(y) + II(x) ■ II(y) 



(4.60) 



The nonlocality arises because the interaction kernel is now a polynomial in |k| rather than 
k. 



V. CONCLUSIONS 

The main, new result of this study was the rigorous derivation of a renormalized effec- 
tive Hamiltonian to order g 2 in the strong coupling constant. This was achieved through 
a similarity transformation in conjunction with a continuous cut-off regularization scheme. 
The use of a sharp cut-off regularization led to an unacceptable non-local Hamiltonian. The 
effective Hamiltonian was determined by introducing an appropriate counterterm Hamilto- 
nian and imposing form invariance on the transformed Hamiltonian. Renormalization not 
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only restored the proper symmetries but also rendered all observables calculated with this 
Hamiltonian insensitive to the cut-off parameter. 

Future formal work will address both the issue of confinement as well as renormalization 
in the combined quark-gluon sector. In particular both perturbative and non-perturbative 
approaches will be examined. Such a treatment is fundamental to understanding the role 
of gluons in the structure of mesons and baryons. Finally, utilizing many-body techniques, 
large scale applications of the effective Hamiltonian focusing upon glueball and hybrid meson 
systems will be reported. 
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APPENDIX: CANONICAL HAMILTONIAN 

In this Appendix we present the canonical Hamiltonian expressed in the Fock represen- 
tation. The free Hamiltonian is simply 



Ha 



d 3 k 



^<t(kK(k) + £ 



d 3 k 



Ek 



b j! (k,s)b(k,s) + d\k,s)d(k,s) 



(Al) 



(2^)3— « — ^ ' (2vr)3 

where E^ = \/k 2 + m 2 and Uk = |k|. The interaction terms will be separated according to 
their term of origin in Eq. (2~T). From the purely gluonic term, H g , we have three-point 
couplings to O(g) 



9 2V2 



l 9 rabc 



j-j- d k n 



fci 



1 1 ., ,,, ^^(2vr) 3 5( 3 )(£k m ) 



(A2) 



x : 



l »(k 1 ) + af(-k 1 ; 



a$(k 2 ) + af(-k 2 ; 



[4(k 3 )+af(-k 3 



where u>x = uj tl , etc. 

For 0(g 2 ) we obtain (normal-ordered) four-point couplings 
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^y(4)_ a s 7r jabc jade 



11 <W 



L x (27r) 3 y v /cJIZJ 2 7^ 3 7J 4 ~ 
oJ(k a ) +af (-k x )] [ S c (k 2 ) +af (-k 2 ) 



(27r)V 3 >CTk r , 



a, d (k,) + af(-k 3 



(A3) 



a|(k 4 ) + af(-k 4 



where a s = g 1 j^. In addition, there are one-body operators which arise in normal ordering 
Eq. QE3|) 

d 3 fcd 3 fc' 1 



H^ = a s nN c J 



(27r) 6 ^fe^fc- 



[2^-A,(k')]e 



-2^ t/ /A 2 



(A4) 



x 



^(k)^(k) + -(<(k)a«(-k)+h.c. 



where N c is the number of colors and a quadratic divergence has been regulated by intro- 
ducing an additional damping exponential. 

The Coulomb term, Eq. (|2.4j) , has a similar structure through 0(g 2 ). To this order we 
may set V • D = V 2 , so that 

p a (x)p a (y) 



H, 



c 



y / d 3 xd 3 y 



+ 0(g 3 ). 



(A5) 



i x - y| 

Writing this in normal order generates two- and one-body operators. We shall exhibit here 
only those terms involving the gluonic component of the color charge density; the purely 
fermionic contributions may be found in Ref. 

The two-body operators involving only gluons are 



o-( 4 >99)_ a s tabctade 
n C 



x: 



h£k) f^iV , V(k 1+ k 2 )(2.)V)(Ek„ 



(kO + af (-kO] [a?(k 2 ) - af (-k 2 



(A6) 



a?(k 3 ) + af(-k 3 )] [a*(k 4 ) 



af (-k 4 



where 



47T 



(A7) 



is the momentum-space Coulomb potential. The one-body operators obtained in normal 
ordering Eq. (A6) are 



H, 



(2) 
C 



a,N, r (k + k>) ( u£±^ D^Me-W 



U k U} k > 



af(k)a|(k) 



(A8) 



a s N c r d 3 kd 3 k' 



(2n) 



V ^k^k' I 



<(k)<(k) + h.c. 
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Again, quadratic divergences have been regulated with exponential cut-offs. 
There are also mixed fermion-gluon two-body operators, given by 



H, 



(4,99) 
C 



ia 



d?k d 3 k' d 3 q d 3 q' 
(2tt) 6 (2tt) 6 



rabcrpa 
J a 



p 



LO n i 



UJr 



1/2 



x 



x 



(27r) 3 5 (3) (-k + k' + q + q') 



a 4 6 (q)+af(-q) 



(A9) 



+ 



t t(k, s)u(k', s')] &£(k, s)^(k', s') + [vt(-k, s)v(-k', s')] tZ a (-k, s)4(-k', s') 
u f (k, s)t;(-k', s')l &£(k, s)4(-k', s') + M(-k, s)«(k / , s')l <U-k, s)6 j3 (k / , s' 



No tadpoles arise from normal ordering this operator since T a has zero trace and the f abc 
are antisymmetric. 

Finally there is the quark-gluon coupling term H qg [Eq. 



H <ig = -9 E 



cflfcdW? TZ, (27r) 3 5 (3 )( _ k + k , + q) ; r a(q) + <t( _ q) 



(A10) 



(2tt)9 

: ([ M f (k, s)a i tt(k , > s')] 6+ (k, s)6 /3 (k / , s') + 
+ [u f (k, s)o i v(-k / , s')l &£(k, s)dt(-k', s') + M(-k, s)a i w(k / , s')l d«(-k, s)6 /3 (k / , s') 



v\— k, s)aif (— k', s') k, s)dt(— k', s') 
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